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Abstract. The stochastic nature of the IMF in young 
stellar clusters implies that clusters of the same mass and 
age do not present the same unique values of their ob- 
served parameters. Instead they follow a distribution. We 
address the study of such distributions, parameterised in 
terms of their confidence limits, in evolutionary synthe- 
sis models. These confidence limits can be understood as 
the inherent uncertainties of the synthesis models. Here 
we concentrate on some parameters such as EW(H/3) in 
emission. For instance, we show that for a cluster where 
10 5 M have been transformed into stars, the dispersion 
of EW(H/3) is about 18% within the 90% confidence levels 
at ages between 3.5 and 5 Myrs. 
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1. Introduction 

Since the development of evolutionary synthesis models 
by Tinsley (1980) several groups have developed models 
trying to reproduce observable quantities in systems in 
which the stellar content is not resolved. These models 
provide a powerful tool to study galactic and extragalactic 
Hii regions. Several types of models can be found in the 
literature. For instance, Leitherer et al. (1996) give an 
extensive and comprehensive review on them. 

The variety of external inputs used in the models (evo- 
lutionary tracks, stellar atmosphere libraries, etc) and the 
treatment of these input parameters, (for example, the in- 
terpolation methods), lead to slightly different model out- 
puts from different authors, although "grosso modo" all 
of them produce similar results. One way to check the ap- 
plicability of these models is to compare the results using 
different inputs. In this direction we highlight the work of 
Bruzual (2000) in which the author presents an extensive 
study on how the external inputs to the models influence 
the output results. 
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These problems can be exacerbated when modelling 
systems with small number of stars given that, in general, 
synthesis models use continuous functions that do not re- 
produce exactly the discontinuous nature of star forma- 
tion, especially in systems with low number of stars. In 
this letter we address this issue by studying the intrinsic 
variations of the evolutionary synthesis models outputs 
due to the deviations of the analytical Initial Mass Func- 
tion (IMF) when modelling Hii regions. 

Here, we study the importance of such fluctuations in 
some typical observables of extragalactic Hii regions. We 
characterise our results showing confidence levels of the 
distribution of a given parameter as a function of age and 
system mass. We focus on the study of the H/3 equiva- 
lent width in emission (EW(H/3), that is extensively used 
as an age indicator of star forming regions. The results 
for the number of ionizing photons, Q(H°), and the ra- 
tio of the blue Wolf-Rayet (WR) bump luminosity to the 
H(3 luminosity are also presented. We differ the study of 
other quantities using different metallicities, IMF slopes, 
star formation regimes and evolutionary tracks to a future 
paper. 

2. Synthesis code and the Initial Mass Function 

The study of the IMF has been broadly covered in the 
astronomical literature. See Scalo (1986) for a compre- 
hensive review. We define the IMF in the following way: 
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where a is the IMF slope, and A is a normalisation factor. 
This function gives us the number of stars in a given mass 
range. The widely used Salpeter IMF corresponds to a = 
2.35 with this definition. The total mass of the system will 
be 
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In an evolutionary synthesis code, the amount of stars 
produced is usually generated binning the IMF (or us- 
ing Montecarlo simulations). The evolution of each bin 
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(or star) is followed along the corresponding evolutionary 
track and the final result is computed from the overall 
evolution of all the stars (or masses) considered. In such a 
procedure, variations in the IMF affect directly the output 
results. The IMF and the synthesised parameters are usu- 
ally normalised to the total mass transformed into stars in 
the synthetic cluster. The comparison with real data is al- 
ways carried out comparing an observed quantity with the 
corresponding one obtained with the code. It is assumed 
that the synthesised quantity is a good approximation to 
the theoretical analytical value and usually no further con- 
siderations are taken into account when comparing data 
and model outputs. 

We have used an updated version of the evolutionary 
synthesis code of Cervino & Mas-Hesse (1994), that uses 
Montecarlo realisations of the IMF to produce star form- 
ing regions of a given mass. We have also used the code 
to generate analytical simulations of clusters of the same 
mass, where analytical means a mass binning approxima- 
tion of the IMF. We have tested that the analytical and 
the Montecarlo outputs show consistent results. The dis- 
persion between them is lower than 1.5% when a large 
number of stars is used (5 x 10 5 stars, i.e. 3 x 10 6 M Q 
transformed into stars in the mass range 2 - 120 M ). 
In this version of the code we use the solar metallicity 
evolutionary tracks of Meynet et al. (1994). The analyti- 
cal results have also been compared with the predictions 
of the code of Leitherer et al. (1999), that uses the same 
evolutionary tracks. The results obtained with both codes 
are in general agreement. 

We have performed 600 Montecarlo realisations of a 
cluster in which 10 3 M are transformed into stars in the 
mass range 2 - 120 M Q ; 400 realisations of a cluster of 10 4 
Mq and 200 of a cluster with 10 5 M Q . The synthesis code 
is the same in all Montecarlo realisations, the only dif- 
ference being the stochasticity of the Montecarlo process 
which reflects the fluctuations of the actual mass distribu- 
tion with respect to the analytical continuous IMF. 

3. Distribution of observed parameters: 
Confidence levels 

In this section we discuss the resulting distributions of the 
observable parameters obtained from the different Monte- 
carlo runs for the masses mentioned above. We note that 
these distribution are in general non-Gaussian. The confi- 
dence levels presented hereafter are computed integrating 
the resulting distributions to contain the 90% (68%) cen- 
tral values of the given parameter. We refer to these val- 
ues of the parameters as the 90% (68%) confidence levels 
of the resulting parameter distribution obtained from the 
Montecarlo realisations. 

Fig. |l|. shows the 90% and 68% (equivalent to ler) con- 
fidence levels for EW(H/3) resulting from the simulations 
in which 10 4 M Q are transformed into stars. This is the 
mass range in which Hn regions with L(Ha) ss 10 39 erg 
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Fig. 1. 90% and 68% Confidence levels of the EW(H/3) 
distribution for the simulations with 10 Mq transformed 
into stars. 

s _1 lie. This value of luminosity of Ha is usually used as 
the boundary between giant and normal Hn regions. 

It can be seen from the figure that the confidence lev- 
els change with time. The largest spread in EW(H/3) is 
obtained between 3 and 9 Myrs (see also fig. |^) when the 
fluctuations in the ionizing flux and the continuum can 
lead to variations of the order of 20% at the la level. The 
dispersion is reduced in the age range between 9 and 13 
Myrs and increases again to an almost constant level for 
older ages. 

The mass dependence of such fluctuations can be seen 
in Fig. H in which the 90% and 68% confidence levels are 
displayed for a fixed age of 4.5 Myr since the onset of the 
burst as a function of the mass transformed into stars. 
The figure shows clearly that the dispersion is not sym- 
metric around the analytical value. The lower confidence 
level shows an almost linear correlation with the logarithm 
of the mass transformed into stars. The upper level does 
not show this linear behavior. An important implication 
regards the different interpretations of observational data 
when compared to models assuming that they provide an 
exact value. As can be seen from the figures, the 90% or 
68% confidence levels of the parameter distributions are 
uncomfortably large, especially for low masses. Due to this 
behaviour of the EW(H/3) evolution, the age determina- 
tion between 3 and 5 Myrs are not well defined at the 
90% confidence level. This situation may be improved us- 
ing more diagnostic observables to estimate the age since 
the onset of the burst. 

Fig. H presents the evolution with age of the 90% 
confidence level of the EW(H/3) distribution for differ- 
ent amount of mass transformed into stars. This figure 
demonstrates the importance of taking into account the 
confidence levels of an observable parameter for a given 
mass range. For instance, for a 10 5 M Q cluster, a value 
of the EW(H/3) equal to 50A, is compatible with any age 
between 3.5 and 5 Myrs within the 90% confidence levels. 
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Fig. 2. 90% and 68% Confidence levels of the EW(H/3) 
distribution for simulations of 4.5 Myrs as a function of 
mass transformed into stars. 
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Fig. 3. Evolution of the 90% confidence level of the 
EW(H/3) distribution for simulations with different 
amount of gas transformed into stars. 

In Fig. |H the evolution of the ionizing flux, Q(H°), is 
displayed against the mass transformed into stars. For a 
10 3 M Q system the fluctuations allow values that differ by 
an order of magnitude for a given age. The 90% confidence 
level values for the distribution of this observable show a 
relatively broad distribution up to 7 Myrs. Comparing this 
figure with the previous one, it is possible to figure out 
that the the scatter in the continuum flux around 4861 A, 
L c (H/3), is the driving parameter producing the time evo- 
lution of the EW(H/3) distribution. We have checked this 
point with the actual values and confirmed it to be the 
case. As a consistency check, we have also corroborated 
that the dispersion in L C (H(3) combined with the disper- 
sion in Q(H°) approximately reproduce the dispersion in 
EW(H/3). 

Fig. ^| shows the evolution of the ratio of the Wolf- 
Rayet bump and H/3 luminosities, L(WRbump)/X(_ff/3), 
presenting the 90% confidence levels of its distribution. We 
would like to point out that the analytical value, close to 
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Fig. 4. Evolution of the 90% confidence level of the Q{H°) 
distribution for simulations with different amount of gas 
transformed into stars. 
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Fig. 5. Evolution of the 90% confidence level of the 
L(WRbump)/L(i?/3) distribution for simulations with dif- 
ferent amount of gas transformed into stars. 



1, is 2.3 times higher that the one obtained in Cervino & 
Mas-Hesse (1994). We have checked this value and confirm 
that it is the result of using different evolutionary tracks 
with enhanced mass loss rates. It is important to note that 
at the 90% confidence level, for a 10 3 M Q cluster, values 
consistent with no detection of WR stars are allowed. The 
distribution at la level is also compatible with such non 
detection. 

This study opens an additional question. If stochas- 
ticity of the stellar mass distribution of low mass clusters 
imply such dispersions in the distribution of observable 
parameters, what would its effect be in a system com- 
posed of several smaller groups? In order to address this 
question, we have simulated multiply-composed 10 5 Mq 
clusters by adding randomly 6 sets of 100 clusters with 
10 3 M and 6 sets of 10 clusters 10 4 M to obtain the ob- 
servables for a variety of ages. We have plotted the results 
for the EW(H/3) and the L (WRbump) / L(H f3) together 
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M Q ; multiplying by 1.5 for the mass range 0.8 to 120 M Q 
and multiplying by 3.75 for the mass range 0.08 to 120 
M . In all the cases a power law IMF with the Salpeter 
exponent has been assumed. 
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Fig. 6. Comparison of the values of EW(H/3) and 
L (WRbump) / L{H j3) obtained for multiply-composed 10 5 
Mq systems (triangles and crosses) with the analytical val- 
ues and 90% confidence level for "single" 10 5 M Q systems, 
(see text) 



with the analytical values and the 90% confidence level of 
the distributions for a 10 5 M cluster in Fig. ^. 

The figure shows that all the simulated points fall 
within the 90% confidence level for the 10 5 Mq simula- 
tion. However, this is not the case for the 68% confidence 
level. In conclusion, the distribution of the observed pa- 
rameters due to the stochasticity of the star formation 
process shrinks towards the analytical values of evolution- 
ary synthesis models as larger masses are considered more 
or less independently of whether these systems are mono- 
lithic unique systems or composed by several coeval sub- 
groups. We would like to stress that our previous simula- 
tions assumed the same age for the subgroups that form 
a system. We have not addressed the consequences of dif- 
ferent ages for the different subgroups. 

As a final note we would like to point out that the total 
masses used here have been computed using stellar masses 
in the range from 2 to 120 M Q , which are the ones relevant 
for the observable parameters discussed. The correspond- 
ing total masses should be corrected by multiplying by a 
factor 1.36 for stellar masses in the range from 1 to 120 



4. Conclusions 

The analysis presented in this letter allow us to draw the 
following conclusions: 

— The intrinsic deviations between analytical and 
stochastic IMF must be taken into account when com- 
paring synthesis models with observational data. The 
inherent uncertainties depend on the amount of gas 
transformed into stars. 

— These deviations are independent of the synthesis code 
and represent a lower limit to the total uncertainties. 
The total deviations will depend on other parameters 
like the star formation history, the input ingredients, 
the numerical approximations, etc. 

— The deviations depend on the evolutionary tracks 
used. 

— The widths of the parameter distributions compared 
to the analytical values are proportional to the mass 
transformed into stars in the stellar cluster or groups 
of clusters. They also depend on the IMF slope and 
the evolutionary status of the cluster. 

— The specific distribution of the deviation from the an- 
alytical value varies from observable to observable. As 
an extreme example, the 68% confidence level of the 
L(WRbump)/L(_ff/3) ratio is compatible with no de- 
tection of WR stars. 
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